In-medium minijet dissipation in Au+Au collisions at y/s^N = 130 and 200 GeV 
studied with charge-independent two-particle number fluctuations and correlations. 
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Medium effects on charged-particle production from minijets are studied using three complemen- 
tary analysis techniques. We find significant angular coUinearity and number correlations on pt 
even at moderate pt < 3 GeV/c. In this pt range abundant particle multiplicities enable preci- 
sion measurements of number correlations of non-identified hadrons for kinematic variables {pt,ri,(f)). 
Methods include (1) direct construction of two-particle correlation functions, (2) inversion of the bin- 
size dependence of non-statistical multiplicity fluctuations and (3) two-dimensional discrete wavelet 
analysis. 

Two-particle correlations on pt exceed expectations from a model of equilibrated events with fluc- 
tuating global temperature. A correlation excess at higher pt is interpreted as final-state remnants 
of initial-state semi-hard collisions. Lower-pt correlations exhibit a saddle structure varying strongly 
with centrality. Variations in the forms and relative strengths of low and high pt correlations with 
increasing centrality suggest transport of semi-hard collision products into the lower pt region as a 
manifestation of in-medium dissipation of minijets. 

Correlations on pt can be associated with angular correlations on {ri,(f)), using analysis methods 
(1), (2) or (3). In particular, wavelet analysis (3) is performed in the {r],(t>) space in bins of pt (< 2 
GeV/c). Observed angular correlation structures include those attributed to quantum correlations 
and elliptic fiow, as well as a localized structure, increasing in amplitude with pt, and presumed 
to originate with minijets. That structure evolves with increasing centrality in a way which also 
suggests dissipation, including an increased correlation length on rj which may be related to the 
influence of a longitudinally expanding medium on minijet fragmentation. 

PACS numbers: 25.75.-q 



This note documents experimental observations of 
medium-modified minijet correlation structures in AuAu 
collisions at RHIC, presented at a poster session of the 
Quark Matter 2004 conference. The data come from 
three distinct analyses which address the same physics 
topic ~ quantitative diagnostic of the strong interaction 
medium created at RHIC. The intrinsic properties of this 
medium in equilibrium are connected with its observable 
response to the excitations , experienced in the course 
of heavy ion collision events due to minijet propagation. 
In this work we analyze the response on the basis of dy- 
namical information contained in fluctuations and corre- 
lations in the number density of non-identified hadrons 
in the space of kinematic variables 77, </> and pt- 

A direct construction of two-particle correlations in 
Pt has been performed on ^/snn ~ 130 GeV data. 
The basic object of this analysis is density of pairs p 
in the two-dimensional space spanned by the pt of the 
two particles. To make statistical errors more uniform 
along the kinematic variable of choice, we replace pt by 
X{pt) = 1 - exp[-(mt - m^)/0.4GeV] € [0,1[. Den- 
sity of sibling pairs (coming from the same event) psib is 
compared with a mixed pair reference Pmix- 

The physics of small momentum differences is known 
to be dominated by quantum statistical and Coulomb ef- 
fects. To focus on the physics of large momentum differ- 
ences (also referred to as large scale) , we eliminate sibling 
and mixed pairs which simultaneously satisfy conditions 
\m-m\ <0.3, |(/>i-</>2| < V6, \pti-Pt2\ <0.15GeV/c 
while having a particle with pt < 0.8 GeV/c. The re- 
sult in form of a differential density ratio is shown in 



Fig. ^ for four centrality classes. Typical magnitudes 
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FIG. 1: Psib/pmix[X{pti),X{pt2)] — 1, pair-density ratios for 
all unidentified charged primary particles for (a) central, (b) 
mid-central, (c) mid-peripheral, (d) peripheral AuAu colli- 
sions at y/SNN = ISOGeV . Note the gray scale differences. 

of Psib/ Pmix ^ 1 Vary with multiplicity as 1/N, but 
there are more subtle shape changes among the pan- 
els of that figure. To quantify those, we form a fit- 
ting model based on a Levy distribution 1 /pt dN/ dpt oc 
[1 -|- /3o("^t — ^ir)]""- Here the exponent n provides a 
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FIG. 2: Left; 2D residuals (data-"soft" fit) for central 
AuAu events projected onto X-e- Shown as a solid curve 
is the "hard" fit based on a hypothesis of a gaussian trans- 
verse rapidity distribution. This hypothesis is consistent with 
pp data. Right: centrality trends in the curvature mea- 
sures. Curves indicate linear trends on mean path length 



parameter responsible for "equilibration" , n —^ oo being 
the Boltzmann limit. It is informative to use the sum 
and difference variables mts = mti + mt2 — 2mjr and 
TOtA = mti — ^12] the ratios in Fig. ^ are concave along 
TOtE and convex along rritA- For the mixed pairs, the as- 
sumed factorization of two Levy distributions results in 
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For the sibling pairs, the fitting model is as in Eq. ^ 
except that n in the first and second term is replaced by, 
respectively, riE and tt-a- Curvatures of PsibI Pmix at the 
origin are /3o(l/nA — l/n)/2 along ttija and (i'^{\/nY. — 
l/n)/2 along mts, and have different signs. 

The strength of the saddle shape is quantified in terms 
of the two curvatures: 
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Thus, the correlation structure is decomposed into sad- 
dle shape (dissipation) and hard component (high X{pt) 
peak). With centrality. the "dissipation" grows with a 
linear trend on the mean path length, as seen in the right 
panel of Fig. |21 

Constructing number correlations directly becomes in- 
efficient when the number of events and multiplicities are 
large, since number of computations scales with the num- 
ber of particles as 0{N'^). To analyze the ^/snn = 200 
GeV data sample, we use fluctuation analyses of compu- 
tational complexity 0{N). 

In the basis of Fourier harmonics, Wiener-Khinchin 
theorem relates autocorrelation with the local fluctua- 
tion power spectrum via Fourier transform. Relations 
between quantities characterizing fluctuations and corre- 
lations can be also obtained in the basis of box functions 
(bins), or in a discrete wavelet basis. 

The inversion of the bin size dependence of the num- 
ber fluctuations is performed to express the results of 




FIG. 3: Reconstructed autocorrelations in particle number in 
AuAu at ^/snn = 200 GeV, for varying centrality. 0a = 
(f>i - <j)2 and rjA =rii - V^- 



the 0{N) fluctuation analysis as a two-point autocor- 
relation (an 0{N'^) quantity). Relation between num- 
ber fluctuation in a ((5?7, 5(j)) bin and net autocorrelation 
AA = Agijj Are/ is 
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where A is autocorrelation. K is kernel which depends 
on the binning used, drj and Scj) are acceptance ranges, 
and 77a, 0a are difference variables. This integral 
equation is solved using standard techniques for inverse 
problems to yield the normalized net autocorrelation 

J^^ ^(?/a,0a), shown in Fig.gl Flow V2 structure 
underlies a same-side minijet peak broadened on tja with 
increasing centrality. 

Application of the DWT power spectrum analysis tech- 
nique in STAR is described in detail in The mea- 
sure of local point-to-point fluctuation is the fluctuation 
power P^{m), constructed out of the DWT expansion co- 
efficients in the (77, 0) space, a^^ij, where A indexes 
the three possible directional sensitivity modes, pseudo- 
rapidity 77, azimuth 0, and diagonal ri(f>. 



P^(m) 
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The basic measure of correlation structure is the so-called 
normalized dynamic texture, (P^rMe ~ Pmix)/ ^mixl^ ^ 
which incorporates the mixed event reference. Here N 
is a number of particles in an event or subevent being 
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FIG. 4: Normalized dynamic texture (see text) as a func- 
tion of pt for central events (0.6 < N/No < 1.1). • - STAR 
data; solid curve - regular HIJING; dashed curve - HIJING 
with jet quenching; □ - peripheral {N/Nq < 0.1) STAR data, 
renormalized according to a "null hypothesis". 
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FIG. 5; Scale de- 
pendence of the dy- 
namic texture (see 
text) in peripheral 
and central events 
for 1.1 < pt < 
1.5 GeV. • - STAR; 
curves are from HI- 
JING. An estimate 
of a systematic error 
due to track merging 
is shown as a hatched 
area. 



analyzed. In case of a onc-dimensional random field X(t) 
such as a stationary time series, the relationship between 
this P{m) and an autocorrelation is: 



P{m)= X{T/2)X{-T/2)W{T,m)dT, (5) 



where t = t2 — ti, and W(t, m) is the kernel function, lo- 
calized and symmetric around 0, which depends only on 
the wavelet. For the Haar wavelet which we use, it is pos- 
itive around and turns negative away from 0, therefore 
one can think of P{m) as a derivative of an autocorrela- 
tion over scale (or r), averaged over an event sample. 

Centrality in the analysis is characterized by the ac- 
cepted number of quality tracks in the TPC, iV, relative 



to A^O: where A^o is such that 99% of minimum bias events 
have N < Nq. 

Fig. ^ shows a difference between dynamic texture 
data for central AuAu events and various expectations. 
In the experiment, we see a change in the pt trend above 
Pt = 0.6 GeV in the 77 mode. Instead of rising with pt 
(as in the peripheral events), the STAR data points be- 
come negative. This contradicts HIJING |^, including 
the "quenching" mode. In HIJING, the rise of the sig- 
nal with Pt is obtained by "turning on" jets. To make 
a comparison between central and peripheral data, by- 
passing the model, we formulate a "null hypothesis" : the 
correlation structure {Ptrue~ Pmix) / Pmix i'^ Au-|-Au col- 
lisions is invariant of centrality. Then, the difference in 
{Ptrue - PLx)/PLx/N between central and peripheral 
events (including the pt trends) is due to the difference 
in 1/A^ (and in dN / dpt) 

Shown in Fig. 0] as □ is the peripheral data from 
STAR, rescaled under an assumption of the "null hy- 
pothesis", taking the difference in dN/ dpt into account. 
In the 77 mode it is visible that the actual and rescaled 
data sets differ in both magnitude and pt trend. This un- 
derscores the difference of correlation structures between 
central and peripheral events and invalidates the "null 
hypothesis" of the 1/A^ scaling. The left panel shows 
that jet-like behavior persists almost unaffected by cen- 
trality in the -qcj) mode. 

We hypothesize that we are observing a modification 
of the minijet structure predominantly in the longitudi- 
nal, 77 direction. Longitudinal expansion of the hot and 
dense medium formed early in the collision makes this 
direction special and is likely to be part of the modifi- 
cation mechanism. If so, we may be observing an effect 
of the longitudinally expanding medium on the minijet 
fragmentation and/or hadronization at "soft" pt. 

Fig. El shows the scale dependence of the correlation 
structure in a pt range where the centrality effect is most 
pronounced. We see that the modification of the struc- 
ture does not imply disappearance of the correlations in 
the central events, even on the fine scales. Thermaliza- 
tion, while affecting the correlations (more than the pop- 
ular model predicts), does not result in a correlation- 
free system in central events, at least as seen by measur- 
ing final state hadrons. 

We conclude that large multiplicities of hadrons in the 
STAR acceptance enable precision studies of the event 
correlation structure. Our observations are consistent 
with an emerging picture: minijets from initial state scat- 
tering arc modified by longitudinally expanding colored 
medium. These measurments of the effect of the medium 
on the parton fragmentation and hadronization provide 
quantitative information about the medium and nonper- 
turbative QCD. 
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